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Abstract 

Pairing plays a central role in nuclear systems. The simplest model for the pairing 
is the constant-pairing Hamiltonian. The aim of the present paper is to include 
the continuum single particle level density in the constant pairing Hamiltonian and 
to make a comparison between the approximate BCS and the exact Richardson 
solutions. The continuum is introduced by using the continuum single particle level 
density. It is shown that the continuum makes an important contribution to the 
pairing parameter even in those case when the continuum is weakly populated. It 
is shown that while the approximate BCS solution depends on the model space the 
exact Richardson solution does not. 
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1 Introduction 

Many-body calculations generally start witli a mean-field wliicli provides tlie 
single-particle representation from which the many-body representation is 
built. The preferred framework for studying nuclear many-body systems is 
the interacting Shell Model (ISM) and its extensions to open systems using 
real [1,2,3] or complex energy [4,5,6] representations. The drawback of the 
ISM is that it becomes cumbersome as the dimension of the model space or 
the valence particles increase. The situation becomes even worse when the 
continuum part of the energy spectrum is included. Even when the pairing 
Hamiltonian includes much less correlations than the ISM it contrives to get 
the most important part of the pairing interaction. The pairing Hamiltonian 
constitutes an important approximation in nuclear systems [7,8,9] and it can 
be considered as a first approximation to a large-scale continuum shell model 
calculation. 
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The pairing Hamiltonian can be solved in an approximate way by using the 
Bardeen-Cooper-Schrieffer (BCS) approximation. In references [10,11,12] the 
BCS solution was studied by using a single particle basis which includes the 
complex energy continuum as well. However the BCS solution has two main 
drawbacks: (i) the many-body wave function does not conserve the number of 
particles and (ii) it does not produces a non-trivial solution for those values 
of the strength which are smaller than a critical value. The constant-pairing 
Hamiltonian can have exact solution [13] if we use the similarity between the 
many-body time-reversed pairs with those of the many-body boson systems. 
The exact solution does conserve the particle number, moreover the solution 
exits for all values of the pairing strength. The eigenf unctions and eigenvalues 
are written in terms of a set of parameters called pair energies. In Rcf. [14] the 
Richardson solution was used in a representation which included the resonant 
part of the continuum. 

In this work an exact solution of the pairing Hamiltonian is given in a repre- 
sentation which includes the continuum in a real energy representation. The 
continuum is included through the continuum single particle level density 
(CSPLD). In order to ehminate the spurious effects of the so called parti- 
cle gas [15,16], the CSPLD relative to a free particle in a box [17] is used. In 
order to asses the limitations of the BCS solution the CSPLD is introduced 
into the approximate BCS solution. 

In Section 2 the pairing Hamiltonian with a continuum basis and the expres- 
sion of the BCS and Richardson equations arc given. The application of the 
methods for the Sn isotopes are presented in Section 3. The conclusion and 
future perspective are summarized in Section 4. 



2 Method 

2. 1 Hamiltonian 

The pairing-model Hamiltonian in a continuum basis reads, 

H = Hsp-G P (1) 
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Here the index /3 = {6,771/3} = {nh^h, 3b,mi3) refers to bound states and 7 = 
{c, m^} = {lc,jc,fn^) to continuum states. The first summation is over the 
valence bound states while the second one is over the partial waves. In practical 
applications upper limits are set for the energy Smax and for the partial waves 
Imax in which continua are taken into account. These parameters determine the 
size of the selected model space. The creation and annihilation a operators 
satisfy the usual anti-commutation relationship with Kroneckcr delta for the 
bound states and Dirac delta for the continuum states. The operator Op — 
(— )-'"~"^''a+_^^ is the time reversed of the a+ operator. The summation u > 
refers to only the positive values of the m^, i.e. projection of the total angular 
momentum. It is assumed that the pairing strength parameter G parameterize 
the total number of particles {A = Acore + ^m/ence) as G = x/A [18]. For a 
given model space the parameter x will be adjusted to reproduce the pairing 
energy taken from theoretical mass table for the isotope of the middle of the 



2.2 Model Space: Bound and Continuum States 

The single particle model space is calculated in a Woods-Saxon (WS) potential 
describing the mean field. For simplicity the energy shift of the levels due 
to the change of the atomic mass number will be ignored in this work, i.e. 
the same single particle energies will be used for all isotopes. The mean-field 
parameters are chosen to reproduce approximately the experimental energies 
of the core plus one nucleon system. The continuum part of the spectrum 
is represented by the continuum single particle level density (CSPLD). The 
CSPLD is defined relative to the free particle density in terms of the phase 
shift 6 [17,15,19,20,21,22,23,24,25,16]. The phase shifts are obtained by solving 
the time independent Schroedinger equation for positive real energy e [26] in 
the WS potential for calculating bound states of the valence particles 



where Sc{e) is uniquely determined by the requirement of its continuity in the 
energy region concerned. 
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2.3 BCS Equations 



The gap parameter A in the continuum basis reads [10], 



A = Ab + Ac 
G 

^ b 

G 7 

Ac = — y de u{e)v{s)g{e) 



(3) 
(4) 



where Ub and u{6) are the usual occupation probabiUties amplitude for bound 
and continuum states respectively in the Biedenharn-Rose phase convention 
[27]. 

The gap and the particle number equations are the following 
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where Ei, and E{e) are the usual quasi-particle energies in the bound and 
continuum states respectively [27] . The CSPLD modifies the pairing strength 
in the gap equation in a effective way and it prevents us from putting non- 
physically large numbers of particles into the continuum. 

The ground state energy Ebcs reads, 



Ebcs = T.i'^Jb + l)vl (^b - -^vl^ + J de v^e) g{e) 
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2-4 Richardson Equations 



A feasible generalization of the Richardson equations to basis with a contin- 
uum part can be found in Ref. [28], 

l-Cj3iM_L_G7,,^ + 2GV — L_ = o (10) 

2 2sb — Ea 2 J 2e — Ea E^ — E^ 

where and E^, are the bound single particle energies in the mean-field and 
the pair energies [29], respectively. The CSPLD is included in the coupled 
equations which determines the pair energies. The number of equations is 
equal to the number of pairs in the valence configuration [29,13]. If we make 
an analytic continuation of Eq. (10) to the complex energy plane it can be 
reduced to the form of ref. [14] . 

The ground state energy written in terms of the pair energies is [29], 

^Rich — '^Ea (11) 
a 

2.5 Correlated energy 

The pair-correlation energy is defined by the difference between the energies 
with and without pairing. The energies with pair correlations in the BCS ap- 
proximation and in the Richardson model are given in Eq. (9) and in Eq.(ll), 
respectively. While the energy without pair correlation in the continuum basis 
reads, 

^ 

The occupation number is the number of particles in the level h below the 
Fermi level A and it is zero above A. The step function ^(A — e) makes the 
integral act only in those cases when the Fermi level is in the continuum. 

The correlation energy reads, 

Eccn-r — Ejyio^i — Eq , (13) 

where model denotes BCS Eq. (9) or Richardson Eq. (11). 
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2.6 Pairing Strength 



The strength parameter x in G = ^, depends on the size of the 'model space' 
and naturally on the 'model solution' too. The model space is define by the 

cut-off energy Emax and the maximum orbital angular momentum Imax, while 
model solution refers to the approximate BCS solution or to the exact Richard- 
son solution. For a given model space the strength parameter % is chosen to 
reproduce the experimental pairing energy calculated from mass table for one 
of the isotopes concerned and it is keep constant for all the other isotopes. 

In order to determine the relationship between the experimental pairing energy 
P and the model-solution pairing energy Pmodei we will use the expression in 
terms of the ground-state energies E [27] and we require P — Pmodei- To be 
specific, let us assume that the proton states form a single close shell and we 
have an even number = 2Npair of neutrons, 

P{N) = 2E{N - 1) - E{N) - E{N - 2) (14) 



For the BCS model-solution Pbcs = A [27] while for the Richardson model- 
solution [13] 



Pmrh — 2£„^r — 



Rich ^'^Pn 

^ pair 



E,,(2N,,,r) (15) 



In order to explain the meaning of the indexes in the single particle energy £p. 
and in the pair energy Ep^ let us consider an example. Let us use the following 

Mb model space: {I/7/2, '2p3/2, Ohg/2, 2pi/2, I/5/2} and Npair = 8, i.e. there 
are = 16 valence particles. We can label the 16 single particle energies in 
the ground state as (£1, 62, ■ ■ ■ , ^le) or as (e^^, Ep.^, . . . , Epg). For the Mb 
model space £p^ — . . . — £p^ = Eif^^^, Ep^ = Ep^ = £^2^3/2' ^pr ~ ^ps ~ 
Then we define the pair energy by the limit lima^Q+Ep. — 2£p.. 

The strength parameter x is fixed for a given isotope for each model-space 
and model-solution. For the BCS x is chosen by the requirement P{N) = A. 
For the Richardson model-solution x is adjusted to reproduce the energy of 
the last pair Re[Ep^^^j2Npair)] = 2£p^^^.^ - P{N). 



3 Application 



The method outlined in secc. 2 will be applied to the Tin isotopes above the 
^^'^Sn core. 
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3.1 Setting the Parameters 



3.1.1 Model Space: Bound and Continuum States 

The mean-field which describes the valence bound states and the continuum 
states is represented by a Woods-Saxon plus spin-orbit potentials with the 
following parameters: Vq = 43.5 MeV, V^o = 13.5 MeV, 

R — ro ^core' '^o — 1-27 fm, a — 0.7 fm. These parameters reproduce approx- 
imately the experimental energies of the single particle states in ^"^"^Sn [30]. 
The bound and the resonant states were calculated by using the computer 
codes [31,32]. The scattering states were calculated by using Ixaru's constant 
perturbation method [33]. 

The bound states for valence orbits are defined by the 2plf0h major shell with 
the following energies I/7/2 = -2.442 MeV, 2^3/2 = -1.395 MeV, O/ig/a = 
-0.923 MeV, 2pi/2 = -0.736 MeV and I/5/2 = -0.148 MeV (in this potential 
the Oii3/2 intruder state is unbound). 

Fig. 1 shows the e„iax dependence of the CSPLD for different l^ax values up to 
^max = 100 MeV. We can easily identify resonances below 20 MeV. They corre- 
spond to the ii3/2, 99/2, ^11/2, ki5/2 and ln/2 partial waves (in increasing energy 
order). We can calculate the widths of the resonances by using the program 
Gamow [31] and get the widths 1(^13/2) = 0.436 x 10"^ MeV, T{gQ/2) = 0.887 
MeV, r(in/2) = 0.382 MeV, 1(^15/2) = 0.251 MeV, T{hj/2) = 1.735 MeV. 
The 530 mesh points in energy used for the calculation of the CSPLD were 
distributed as follows: 50 mesh points were taken between zero and the first 
resonance. Then 10 mesh points were taken for each resonance. Above 20 MeV 
one mesh point /MeV was used. 




10 20 30 40 50 60 70 80 90 100 
e [MeV] 



Fig. 1. CSPLD up to Emax = 100 MeV for different Imax values. 
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An increase of the Imax value makes the potential less important and 5c{s) — > 
as Ic ^ oo for a fixed energy and as e — )■ oo for a fixed angular momentum 
[26]. This implies that increasing l^ax will affect mostly the tail of g{s) and will 
not affect the low energy spectrum in which we are interested in. Therefore it 
will be enough to take Imax — 20. 

We will consider two model spaces: (i) Mf,: formed by only the valence bound 
states, and (ii) Mc'. formed by all partial waves up to 20 with a cutoff energy 
at 100 MeV. 



3.1.2 Pairing Strength 

We choose the nucleus ^'^^Sn, with = 16 as the reference nucleus to fit the 
strength parameter x- For this nucleus there is no experimental data available 
therefore theoretical masses were used instead [34,35]. The pairing energy 
obtained was P(16) = 0.985 MeV. Solving the BCS equation for the model 
space Mb gave xf*"*^ = 23.66 MeV, while for the model space Mc, Xc'~^^ — 
14.67 MeV. Solving the Richardson equations we obtained x^^^'^ — 18-4 MeV, 
and xf*^'^ = 12.25 MeV. We observed a much smaller pairing strength value 
in the Richardson model compared to the BCS model. 

A complete representation should include the continuum part of the energy 
spectrum. The absence of the continuum in the model space Mf, shows up in 
a larger value of the pairing strength x- 



3.2 Results 

After fixing the model space and the paring strength we can calculate the 
quantities introduced in the section 2. 

3.2.1 BCS model solution 

The gap parameter A was calculated in both Mj, and M^ model spaces for the 
tin isotopes between ^^^Sn and ^^^Sn. Their values are shown in fig. 2 

In the Mc model space the pairing gap is smaller if N is below A" = 16 and 
it is larger for A" > 16. As a consequence the isotopes being close to the 
drip line have a finite gap value, while with the Mb model space they go to 
zero rapidly. In order to check the independence of the previous conclusion on 
the cutoff energy the gap parameters were calculated for a model space with 
Ejnax — 30 MeV (x = 15.58 MeV) also as shown in the same figure. 
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Fig. 2. Gap parameter for Tin isotopes in the bound model space (Mb), and in two 
continuum model spaces: Emax = 100 MeV (Mg) and Smax = 30 MeV. 

Fig. 3 shows the contributions of A;, Eq. (3) and Ac Eq. (4) to the constant gap 
A in the model space Mc- It shows that the continuum part of the spectrum 
makes an important contribution to the pairing gap. 
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Fig. 3. Relative contribution to the gap parameter for Tin isotopes in the model 
space Mc- 



Because the continuum has an important contribution to the pairing, we can be 
interested in the occupation of the states in the continuum. Table 1 shows the 

contribution of A^^ Eq. (7) and N^, Eq. (8) to the particle number A^. It shows 
that the "occupation" of continuum states does not exceed 4%. The same 
figure was found when we took a different cutoff energy. This a consequence 
of using the CSPLD relative to the free particle density. 
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Finally we calculated the correlation energy Eq. (13) for the two model spaces 
Ml, and Mc as it shown in figure 4. It is found that above ^^^Sn {N = 14) the 
continuum increases the pairing correlation. 

3.2.2 Richardson model solution 

The Richardson model gives the exact solution to the constant pairing Hamil- 
tonian. The pairing correlated energy was calculated using the Richardson 
model solution in the model space without continuum Mf, and the model 
space with continuum Mc- Figure 5 shows the result. The solutions in the two 
models do not depart much from each other. 

3.2.3 Comparison of the solutions in the BCS and Richardson models 

The BCS model solution is an approximate solution of the pairing Hamilto- 
nian which does not conserve the particle number, while the solution in the 
Richardson model is the exact for the constant pairing Hamiltonian. In spite of 

Table 1 

Contribution from the continuum (TVc) and bound (iV;,) particle number TV = 
Nc + Nb for the model space Mc. 



N 


Nb 




2 


1.917 


0, 


.083 


4 


3.848 


0, 


.152 


6 


5.787 


0, 


.213 


8 


7.718 


0, 


.282 


10 


9.619 


0, 


.381 


12 


11.504 


0, 


.495 


14 


13.391 


0, 


.609 


16 


15.285 


0, 


.715 


18 


17.189 


0, 


.811 


20 


19.102 


0, 


.898 


22 


21.026 


0, 


.974 


24 


22.955 


1, 


.045 


26 


24.886 


1, 


.114 


28 


26.808 


1, 


.192 


30 


28.498 


1, 


.502 
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Fig. 4. Pairing correlated energy in the BCS model solution for the M;, and 
model spaces. 
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Fig. 5. Pairing correlated energy in the Richardson model solution for the and 
Mc model spaces. 

this huge difference between the two models the ground state energies for the 
Tin isotopes ^^^Sn-^'^^^Sn fig. 6, shows small differences between the results of 
the two approaches for the isotopes there. A small difference appears between 
the BCS and Richardson solutions after the single particle ground state is 
filled completely. A small difference between the two model spaces in the BCS 
solution can be observed from N — 12 on. The variational BCS solution gives 
a smaller energy than the exact Richardson solution. Since the BCS is a varia- 
tional solution its energy should be greater than the exact one. This apparent 
contradiction is resolved if we remember that the solutions were calculated 
with different value of the pairing strength. The clue of the good agreement 
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is that the binding energy is an observable and it should not depend neither 
on the model solution nor on the model space. The differences between the 
two model solutions and between the model spaces were absorbed by a single 
parameter, i.e. the strength parameter Even when this last statement is a 
trivial affirmation from a theoretical point of view it is not trivial in practical 
calculations. 




N 

Fig. 6. Comparison of the ground state energy for the Tin isotopes in the approx- 
imated BCS solution and the exact Richardson solution for the two model spaces 
Mh and Mc. 



In the BCS model calculation the pairing-correlated energy shows an enhance- 
ment of the pair correlation as the Fermi level A approaches the threshold. The 
Richardson model solution does not show a similar trend if we change from 
the Mb model space to the Mc model space. This shows that when the contin- 
uum is included the BCS approach artificially enhances the pair correlations 
for a Fermi level being less or around 1 MeV (for N — 16 the Fermi level is 
A = -1.12 MeV). 

The inclusion of the continuum in the BCS solution predicts that the pairing 
energy will not be zero at the end of the shell as it happens in a bound 
representation. The Richardson solution confirms this result with a larger 
value than the one predicted by the BCS. 

By comparing the paring-correlated energies of the two models we found that 
the BCS overestimates the pair correlation for every isotopes of the chain 
(except the last one) by an average amount of 278 keV. The largest difference 
is 848 keV at the middle of the isotope chain. 
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4 Conclusion 



In this paper we calculated the many-body constant-pairing model Hamil- 
tonian in the framework of the approximate BCS solution and in the exact 
Richardson solution. A representation which included the continuum part of 
the spectrum using the level density was used. The level density was defined 
with respect to the box normalization in order to subtract from the total den- 
sity the contribution coming from the Fermi gas. We have found that even 
though the continuum makes an important contribution to the pairing gap, 
the population of the continuum is weak and its value is independent of the 
cutoff energy. 

At the scale of the binding energy we have shown that the ground state energy 
is almost independent of the models considered and it is also indcpcnt on the 
size of the model space. But at the scale of the pair correlated energy the 
BCS solution does depend on the size of the model space. This dependence is 
enhanced for Fermi energies being around 1 MeV below the threshold where 
the continuum starts. The correlated energy is overestimated by an average 
amount of 300 kcV in the approximation which does not conserve the particle 
number. We demonstrated that with the exact solution the results are almost 
independent of the size of the model space. This poses the question why to 
include the continuum part of the spectrum into the basis. This question can 
not be answered with the observables calculated in this paper. It is worth to 
make an observation at this point, we should include the continuum part of the 
spectrum into the basis because doing so the effective parameter dependence 
(the paring strength dependence in the constant pairing case) is reduced. 

The future perspective includes the study of the spectra of the even Carbon 
isotopes in the Richardson model solution with continuum single particle level 
density and the comparison of the real energy representation with the complex 
energy representation of the continuum spectrum. 

This work has been partially supported by the National Council of Research 
PIP-77 (CONICET, Argentina). 
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